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or  TAPERED  CANTILEVER  BEAMS  (U) 


6.V.  Price 


1 . Introduction 

The  Defence  Research  Establishment  Suffield  (ORES),  in  support 
of  Canadian  Forces  (Maritime)  requirements,  is  conducting  a series  of 
tests  to  determine  the  ability  of  certain  antenna  designs  to  withstand 
blast  overpressures  of  various  intensities.  Two  of  the  antennas  to  be 
evaluated  are  a 35  ft  fibreglass  Whip  Antenna  and  a 23  ft  UHF  Polemast 
Antenna.  The  objective  of  this  study  is  to  develop  a numerical  procedure 
to  predict  the  time  response  of  antennas  of  the  above  type  when  subjected 
to  a transient  air  blast  load. 

The  numerical  procedure  begins  with  the  Bernoulli -Euler  equation 
for  a tapered  cantilever  beam  subjected  to  a transient  distributed  force. 

This  equation  is  based  on  linear  elastic  theory  and  assumes  small  beam 
deflections  [1].  The  normal  modes  and  natural  frequencies  of  the  beam  are 
determined  by  solving  the  difference  equations  for  free  vibration  using 
successive  relaxation,  Rayleigh  quotient  and  Gram-Schtnidt  orthogonal ization 
numerical  techniques  [2,3].  The  forced  vibration  solution  is  determined 
using  normal  mode  coordinates  and  Laplace  transforms  [4].  In  this  calculation. 
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the  transient  air  blast  load  is  computed  using  the  classical  gas  dynamics 
equations  describing  air  shock  waves  [5,6,7]  and  recent  empirical  drag 
loading  relations  developed  at  DRES  [8,9]. 

The  Bernoulli -Euler  equation  for  a vibrating  beam  of  length  L may 
be  written  in  the  form  [1] 

pA  + -3L[ei-^]=  w,  OsxsL,  (1) 

3c2  \ 3x2  j 

where  x is  distance  from  the  base  of  the  beam,  t is  time,  y(x,t)  is  the  beam 
deflection  in  the  transverse  direction,  p is  the  beam  density,  A(x)  is  the 
beam  cross-sectional  area,  E is  Young’s  modulus  for  the  beam  material,  l(x) 
is  the  second  .romcnt  of  U-e  cross-scctional  area  with  respect  to  the  neutral 
axis,  and  w(x,t)  is  the  load  per  unit  length  on  the  beam  ip  the  transverse 
direction.  The  discussion  which  follows  considers  a numerical  solution  to 
the  above  equation  at  N prescribed  points 

x^  = i AX,  AX  = LN"^,  1 =5  1,  ....  N.  (2) 

subject  to  specified  initial  conditions  and  boundary  conditions. 

2 . Free  Vibration  of  a Tapered  Cantilever  Beam 

The  normal  modes  and  natural  frequencies  are  obtained  by  solving 
equation  (1)  with  w=0.  Assuming  a sinusoidal  time  response,  the  free 
vibration  solution  may  be  written  in  the  form 

y(x,t)  = Y^'"^(x)e^'^r^  , (3) 
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/ y'  \ 

where  w is  the  r'th  circular  frequency  of  the  time  response,  Y'  ^(x)  is 
r 

the  normal  mode  corresponding  to  and  i is  >^.  With  the  results  (2) 
and  (3),  equation  (1)  with  w=0  reduces  to  the  form 


pA 


. Ei, -Siy*''), 

dx2  \ ’ dx^  M 


i 'i 


(4) 


at  point  . Expanding  the  derivative  in  (4)  using  central  finite 
differences,  the  difference  equation  for  free  vibration  becomes 


^i-i  A'-l  ^ " ('i-i " ^li  ^ ^1+1  - 


-2(1,  + I...,)YSr]  + I.O  =*  0,  i = 1.  2,...,N. 


(r) 


(5) 


•T I 111 


1*1  I ' k 


The  four  unknowns  Y_^  , Yq  , and  Y|^^^  are  eliminated  from  this  system  of 
equations  through  the  application  of  boundary  conditions.  Several 
possibilities  are  considered  below 
(a)  clamped  at  x = 0,  zero  displacement  and  slope: 

y(0)  - 0 , Y^"'^  = 0 ; ] 


9X 


= 0 


y{r)  _ yCr) 
^-1  ■ ^1  • 


(6) 


(b)  free  at  x = L,  zero  moment  and  shear: 

£x.' 


9x2 1 


3x3 


= 0 


= 0 


’n+1 


M . 


’N+2  N-1  ^N-2 
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(c;  pini-ied  at  x = 0 and  x = I Ax,  zero  displacement  and  moment  at  x = 0, 
and  zero  displacement  at  x = I ax: 


3x2 


^(0)  = ^ 

y{iAx) 


= 0 . 


^ . y(r)  _ (r) 


y{'‘1  = 0 


I (8) 


With  a suitable  combination  of  the  boundary  conditions  (6)  to  (U), 
equation  (5)  becomes  an  eigenvalue  problem  for  an  n-rowed  square  matrix. 

In  matrix  form,  the  eigenvalue  problem  may  be  written  in  the  form 


0)^  e Y 


ir) 


= a ¥ 


Ar) 


(9) 


where  a is  a f ive-diagorta'  matrix  which  depends  on  the  coefficients 

i=l,...N,  and  b is  a diagonal  matrix  for  which  thei'th  diagonal  element 
is  the  value  pA^ax'*L  '^.  If  the  frequencies  are  all  distinct,  it  may 
be  readily  verified  that  the  normal  modes  are  orthogonal  with  respect  to 
the  matrix  b,  according  to 


y{r}  g y(s] 


= <5  M 

ir<:  r 


0 , r=ts  , 

1 , r=s  . 


(10) 


A numerical  solution  procedure  for  the  eigenvalue  problem  (9) 
will  now  be  considered.  The  procedure  begins  with  an  init'.al  gues  for  the 
smallest  circular  frequency  u^.  The  corresponding  normal  mode  is  obtained 
by  successive  relaxation  [2,3]  according  to  the  iteration  equation 
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Vj''  = Uf  Oij  vj 


1 

(n+1  guess) 


j=l 

(jfi) 


(1) 
j 

(n  guess) 


(11) 


This  iteration  is  inherently  unstable  since  the  system  of  equations  is 
homogeneous.  However,  convergence  to  a stationary  ^ is  achieved  by 
normalizing  after  each  complete  repetition  of  the  iteration.  An 
improved  estimate  for  is  then  obtained  by  the  Rayleigh  quotient  [2,3] 
technique. 


hi 


N 


“1 


T.  M)  I!  a yC) 

1=1  M ^'i1  j=l  ij 

^ v(1) 

L 'i  h. 

i=l 


(12) 


The  complete  iteration  is  repeated  until  a stationary  value  for  the 
circular  frequency  w-i  is  obtained,  at  which  point  the  corresponding 
nt.'-wal  mode,  in  normalized  form,  will  be  Y^^^  . 

(r) 


Hiqher  natural  freauencies  w and  normal  modes  Y' 

r 


are  obtained 

using  the  above  procedure  with  but  one  modification:  between  repetitions 


of  the  successive  relaxation  iteration  (equation  11),  the  latest  guess 
(r)* 

Yj  is  modified  to  extract  any  cotiiponent  in  the  direction  of  previously 

tsl 

evaluated  normal  modes  Y'  s<r.  This  operation  is  termed  Gram- Schmidt 


orthogonal ization  [2,3],  and  is  based  on  the  orthogonality  property  (10). 

r-1  ■ 

y('’)  , Y(r)’‘  ^ ^ ^(s) 

(new  guess) 


^y('^)*, 


(13) 
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In  this  equation,  the  inner  product,  denoted  <u,  v>  for  two  arbitrary 
vectors  u and  v,  is  weighted  with  respect  to  the  matrix  e according  to 

<u,  v>  = u^6v  , (14) 

where  superscript!  denotes  the  transpose  of  the  indicated  vector.  A 
computer  program  which  determines  the  natural  frequencies  and  normal  modes 
according  to  the  techniques  described  above  is  given  in  the  Appendix. 


3"  Forced  Vibration  of  a Tapered  Cantilever  Beam 

The  forced  vibration  solution  begins  with  an  expansion  of  the 
displacement  ve(,tor  y in  terms  of  normal  mode  coordinates  n,  according  to 


y 


N 


z 

r=l 


or 


y = >!>  n , 


where  * is  a matrix  formeu  by  the  normal  modes. 


\ 


/ (15) 


/ 


4 


(1) 

y(n) 

... 

1 

(1) 

HI 

y(n) 

1 HI 

(16) 


Using  the  general  techniques  outlined  in  [4],  substituting  (15)  into  (1) 
and  pre-multiplying  the  resulting  equation  by  4^,  the  governing  equation 
becomes 
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T ••  T T 

4“  + 4 a4n  = 4 W » 


(17) 


where  " denotes  the  time  derivative  dVdt^.  The  r'th  equation  in  the 
system  of  equations  (17)  is 

6V(''>  -n . y a y''">  n,  - w . 

Since  the  normal  modes  satisfy  equations  (9)  and  (10),  the  above  result 
simplifies  to  the  form 


(18) 


,(r)' 


••• 


w(t) 


(19) 


m: 


The  normal  mode  coordinates  nj.vr=l,...N  are  obtained  by  solving 
equation  (19)  using  Laplace  transforms  and  the  convolution  theorem.  With 
specified  initial  conditions  in  the  form 


y(x,0)  = f(x)  , 

y(x,0)  = g(x)  . 

The  solution  to  equation  (19)  becomes 


(20) 


• 1 


(rl^ 

n^(t)  - Y'  ' (pf  cosw^t  + eg  sinu^t 


sinu^(t-T)dT) 

n 


(21) 


The  parameter  appearing  in  this  equation  is  evaluated  according  to 
equation  (10). 


(22) 
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Combining  the  results  (21)  ond  (15),  the  beam  displacement  at  position 


becomes 


y^(t) 


N 

1 N \ 

Z h;' 

r-l  ' y 

COS(i3„t 

r 

L 

/ ^ 

+ \'sin..^t  I pA^gj 


(23) 


The  time  integral  in  (23)  may  be  readily  evaluated  by  numerical  integration 
using  a time  step  At  of  specified  duration. 

A computer  program  which  determines  the  forced  vibration  solution 
according  to  the  techniques  described  above  is  given  in  the  Appendix. 


i 

! 


[, 

f 

\ 


^ • Transient  Drag  on  Circular  Cylinders  from  Air  Blast  Waves 

The  accuracy  of  the  forced  vibration  solution  is  strongly  dependent 
on  the  procedure  used  to  calculate  the  air  blast  load  on  the  tapered  canti- 
lever beam.  The  transient  air  blast  load  on  circular  cylinders  has  been 
investigated  extensively  at  ORES  over  the  period  1968  to  1975.  The  discussion 
which  follows  is  based  on  the  more  recent  DRES  reviews  of  empirical  relations 
for  the  prediction  of  drag  loading  [8,9]  and  on  the  classical  gas  dynamics 
equations  describing  air  shock  waves  [5,6,7]. 
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The  air  blast  load  arising  from  the  Interaction  of  a blast  wave 
with  a circular  cylinder  may  be  expressed  in  the  form 

F = P*A  + Cq  qj  A . (25) 


where  A is  the  projected  area  of  the  cylinder,  P*  is  an  empirical  pressure 
parameter  which  accounts  for  the  load  on  the  cylinder  over  the  short  period 
of  time  during  which  the  blast  wave  diffracts  around  and  engulfs  the  cylinder, 
Cq  is  an  empirical  aerodynamic  drag  coefficient,  and  Oj  is  the  Impact  pressure 
(instantaneous  difference  between  the  local  stagnation  and  static  pressures). 

The  following  empirical  relations  have  been  determined  at  DRES  for 
the  diffraction  pressure  P*  and  the  drag  coefficient  Cq  [8,9]: 


P*  - { 


-1  1.13 

1.6  tj  (PQtj) 


0<t<5t 


T ’ 


0 


t>5t. 


(26) 


M>0.48,  Re  >3x10^ 
M<0.48,  Re  >3x10^ 
M<0.48,  Re  <3x10^ 


(27) 


0.7 

0.6 

I 

In  the  above  equations,  p^  is  the  peak  overpressure  of  the  blast  wave,  tj  is 
the  transit  time  required  by  the  blast  wave  to  completely  engulf  the 
cylinder,  M is  the  instantaneous  Mach  number  of  the  flow  incident  on  the 
cylinder,  and  Re  is  the  instantaneous  Reynolds  number  of  the  flow  incident 
on  the  cylinder. 

To  complete  the  transient  drag  calculation,  it  remains  to  evaluate 
the  basic  fluid  properties  of  the  flow  behind  the  blast  wave.  The  parameters 
to  evaluate  include  q^,  ty,MandRe.  The  evaluation  of  these  parameters  requires 
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5 preliminary  evaluation  of  the  temperature  T,  speed  of  sound  c,  static 
p.-essure  P,  dynamic  pressure  q.  fluid  particle  velocity  u,  and  fluid  viscosity 
0 both  at  the  cylinder  and  immediately  behind  the  shock  front.  To  simplify 
the  presentation,  the  following  notation  will  be  adopted.  Variables  which 
apply  limedlately  In  front  of  the  shock  front  will  be  indicated  by  a subscript 
"A",  variables  which  apply  inmediately  behind  the  shock  front  will  be  Indicated 
with  a subscript  "B",  and  variables  which  apply  instantaneously  at  the  cylinder 
Itself  will  have  no  explicit  subscript. 


Shock  Front 


Cylinder 


The  pressure-time  history  of  the  ideal  olast  wave  is  most  frequently 
described  in  terms  of  the  modified  Friedlander  equation  [7], 

P(t)  = Po(l-t/tj)e’'"^/^d  , Oststjj  . (28) 

where  p(t)  is  the  blast  wave  overpressure,  t is  the  time  measured  from  time 
of  arrival  of  the  blast  wave,  is  the  peak  overpressure,  t^  is  the  positive 
phase  duration,  and  k Is  the  exponential  decay  constant.  The  three  parameters 
Pq,  tjj  and  < in  this  equation  make  it  possible  to  match  the  following  three 
experimental  characteristics  of  the  real  blast  wave:  p^,  t^  and  1^,  where 
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is  the  positive  phase  impulse  defined  by 


f P(t) 
-'0 


The  region  immediately  in  front  of  the  shock  front  is  assumed  to 
be  in  an  undisturbed  state  (u^=0)  at  a known  pressure  and  temperature  T^. 
Conditions  immediately  behind  the  shock  front  are  determined  from  the 
Rankine-Hugoniot  and  other  standard  gas  dynamics  equations  describing  shock 
waves  in  air.  Assuming  a specific  heat  ratio  y = 1-4,  and  adopting  the 
notation  5 to  represent  the  pressure  ratio  Pg/P^  across  the  shock  front,  the 
fluid  properties  immediately  behind  the  shock  front  may  be  written  in  the  form 


Pb  = Pa^Po 


Tg  = C (6  + t)(l  + 6c)-i  , I (30) 

Us  = (1.4  R T^)i(l  + 6p^/(7P^))^  , 

whore  R is  the  gas  constant  for  dry  air,  and  Ug  1s  the  s peed  of  the  shock  front. 
At  the  cylinder  itself,  conditions  are  determined  by  assuming  isentropic  flow  in 
the  region  behind  the  shock  front. 

P = P;\  + P(t)  , 

T = Tg(P/Pg)-286  , 

C = (1.4  RT)^  , 


q = 2.5  p^2(pq  + 7P^)-i(p(t)/Po)2  . 

P “ P(RT)-i  , 


u = (2q/p)  * . 
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With  the  basic  fluid  properties  known  behind  the  shock  front  and  at 
the  cylinder,  it  is  possible  to  evaluate  the  four  parameters  qj,  tj,  M and  Re 
appearing  in  the  drag  equations  (25)  to  (27). 

M = u c"^  , 

Re  = u D 

(32) 

t-j.  = D u^'i  . 

qj  = q + .175  PM**  + .0175  PM^  , 

where  the  kinematic  viscosity  v is  a known  function  of  T and  P (refer  to  [10] 

and  [11]  for  considerations  in  this  regard). 

A computer  program  which  determines  the  air  blast  load  on  circular 
cylinders  according  to  the  techniques  described  above  is  given  in  the  Appendix 

5 . Numerical  Predictions  and  Comparison  with  Experiments 

The  numerical  simulation  model  developed  above  was  used  to  predict 
the  time  response  of  a 35  ft  fibreglass  Whip  Antenna  and  a 23  ft  UHF  Polemast 
Antenna  when  subjected  to  a nominal  7.0  psi  peak  overpressure  blast  loading. 
The  two  antennas  were  subsequently  exposed  to  air  blast  loading  in  Event  Dice 
Throw,  a 620  ton  AN/FO  free-field  blast  trial  conducted  by  the  United  States 

Defense  Nuclear  Agency  at  the  lihite  Sands  Missile  Range  in  New  Mexico  on 

October  6,  1976,  and  strain  data  were  recorded  for  five  pairs  of  strain  gauges 
measuring  bending  strain  on  each  antenna.  The  theoretical  predictions  and 
experimental  results  for  the  Whip  and  Polemast  Antennas  have  been  examined  in 
detail  in  references  [12]  and  [13].  The  brief  discussion  which  follows 
considers  only  representative  items  from  the  detailed  reports. 

The  s*-ructuresof  the  antennas  were  represented  in  the  computer  model 
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in  such  a way  as  to  simulate  the  mass  and  projected  (normal  to  blast  direction) 
cross-sectional  area  profiles  of  the  prototype  antennas.  Photographs  of  the 
two  antennas  installed  at  the  nominal  7.0  psi  peak  overpressure  location  for 
Event  Dice  Throw  are  shown  in  Figures  1 and  2.  The  physical  features  which 
describe  the  corresponding  computer  simulation  of  the  antennas  are  outlined  in 
Tables  1 and  2.  It  should  be  noted  that  the  computer  simulations  of  the 
antennas  agree  with  the  actual  structures  in  the  following  critical  areas: 
weight  distribution,  total  weight,  projected  (normal  to  blast  direction) 
cross-sectional  area  distribution,  and  total  projected  cross-sectional  area. 

Prior  to  the  blast  trial,  a static  load  was  applied  near  the  top  of 
each  antenna  using  an  anchored  steel  cable  at  a pull  angle  of  30°  to  the 
horizontal.  The  load  was  subsequently  released  electrically  ("Twang  Test", 
see  [12,13])  and  the  strain  data  for  free  vibration  were  recorded.  A Fourier 
analysis  was  performed  for  tne  experimental  strain  data  to  determine  the 
natural  frequencies  of  the  antennas.  A comparison  of  theoretical  (numerical 
simulation)  and  experimental  ("Twang  Test")  natural  frequencies  for  the  two 
antennas  is  presented  in  Table  3.  It  is  apparent  from  this  comparison  that 
the  predicted  frequencies  are  in  good  agreement  with  the  values  obtained 
experimentally. 

A comparison  of  theoretical  (numerical  simulation)  and  experimental 
(blast  trial)  bending  strain  histories  is  presented  in  Figures  3 and  4.  The 
theoretical  predictions  were  generated  using  a Friedlander  overpressure  v.ave 
(Figure  5)  which  corresponds  to  the  Defense  Nuclear  Agency  (DNA)  pre-trial 
predictions  for  peak  overpressure  (7.0  psi),  positive  duration  (242  msec)  and 
positive  phase  Impulse  (600  psi-msec).  The  free-field  overpressures  at  the 
base  of  the  antennas  were  measured  with  four  strain-type  pressure  transducers. 
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The  Friedlander  overpressure  wave  which  corresponds  to  the  averagf  experiniental 
peak  overpressure  (6.6  psi),  positive  duration  (251  msec),  and  positive  phase 
impulse  (705  psi-msec)  is  shown  in  Figure  5 for  comparison  purposes. 

It  is  apparent  from  Figure  3 that  the  Whip  Antenna  bending  strain 
predictions  are  generally  much  larger  than  the  corresponding  experimental 
strains  (the  differences  are  within  experimental  error  bands,  considering 
pressure  transducer  errors,  drag  coefficient  errors,  material  uncertainties, 
etc.).  The  average  ratio  of  peak  theoretical  to  experimental  bending  strain 
from  all  five  Whip  Antenna  gauge  pairs  is  1.65  (Table  4).  This  may  be 
compared  to  the  results  in  Figure  4 for  the  Polemast  Antenna  in  which  there 


is  excellent  agreement  between  the  theoretical  and  e.xperimental  strains.  The 


average  ratio  of  peak  theoretical  to  experimental  bending  strain  from  all 
five  Polemast  Antenna  gauge  pairs  is  1.19  (Table  4).  A considerably  more 
detailed  evaluation  of  the  above  two  antenna  experiments  may  be  found  in 
references  [12]  and  [13]. 


I u;>  iv^iia 

A numerical  procedure  was  developed  to  predict  the  elastic  response 
of  variable  cross-section  cantilever  beams  when  subjected  to  a transient  air 
blast  load.  The  numerical  model  was  used  to  predict  the  time  response  of  a 
35  ft  fibreglass  Whip  Antenna  and  a 23  ft  UHF  Polemast  Antenna  when  subjected 
to  a nominal  7.0  psi  peak  overpressure  blast  loading.  The  computed  natural 
frequencies  and  transient  strains  were  in  reasonable  agreement  with  value:- 
obtained  experimentally  in  Event  Dice  Throw.  A more  detailed  evaluation  of 
the  above  two  antenna  experiments  may  be  found  in  references  [12]  and  [13]. 
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X (in) 

ID^in) 

0D*(1n) 

OD^(in) 

0.00** 

4.400 

6.50'’' 

6.5003 

41.85 

4.288 

5.942 

5.942 

83.70 

4.002 

4.878 

4.878 

125.55 

3.665 

4.459 

4.459 

167.40 

3.403 

4.154 

4.154 

209.25 

3.056 

3.935 

3.935 

251.10 

2.707 

3.241 

3,241 

292,95 

2.298 

2.63^ 

2.638 

334.80 

1.957 

2.257 

2.257 

376.65 

1.695 

2.030 

2.030 

418.504 

1.500 

1.900^ 

1.9003 

This  profile  establishes  the  mass  distribution, 

2 

This  profile  establishes  the  pro.lected  (normal  to  blast  direction) 
cross-sectional  area  distribution. 

3 

Extrapolated  value  based  on  data  supplied  by  the  manufacturer. 

Boundary  conditions;  clamped  at  x=C  in,  free  at  x=41S.50  in. 

E = 3.9xl06psi 

P >=  0.002298  slugs/in3 

AX  = 41 .85  in 

N = 10 

L = 418.5  in  (34.88  ft) 

= 12.58  psi 

T^  = 54.0°F 

Po  = 7.0  psi 

tjj  = 242  msec 

Ij  = 600  psi -msec  (k=1.137,  computed) 

At  = 1.00  msec 

The  time  response  is  formed  using  only  the  lowest  3 natural  frequencies 
and  corresponding  normal  modes. 

table  1 : Numerical  model  simulation  of  a 35  ft  fibreglass  Whip  Antenna. 

All  symbols  are  defined  in  the  text.  The  antenna  outside 
diameter  (OD)  and  inside  diameter  (ID)  data  were  interpolated 
from  data  supplied  by  the  manufacturer.  The  air  blast  predictions 
were  provided  by  the  United  States  Defense  Nuclear  Agency. 
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X (in) 

ID^in) 

OD^in) 

OD^in) 

03 

8.978 

9.500 

9.500 

36® 

8.978 

9.500 

9.500 

72 

8.978 

9.500 

9.500 

108 

7.950 

9.500 

17,220 

144 

8.978 

9.500 

9.500 

180 

8.600 

9.500 

13.360 

216 

8.600 

9.500 

13.360 

252 

8.290 

9.500 

10.080 

288® 

8.290 

9.500 

10.080 

This  profile  is  calculated  to  establish  the  correct  mass  distribution, 
assuming  a fixed  OD  equal  to  that  of  the  seamless  extruded  aluminum 
tubing  which  constitutes  the  primary  structural  portion  of  the  antenna. 

2 

This  profile  is  calculated  to  establish  the  correct  projected  (normal 
to  blast  direction)  cross-sectional  area  distribution. 

3 

Boundary  conditions:  pin  at  x=0  in,  pin  at  x»36  in,  free  at  x=286  in. 

E = 10x10®  psi 

p = 0.003044  slugs/in^ 

Ax  = 36.0  in 

N = 8 

L = 288,0  in  (24.0  ft) 

Pp  = 12.58  psi 

= 54.0’F 

,'o  = 7.0  psi 

= 242  msec 

Ij  - 600  psi-msec  (k=1.137,  computed) 

At  = 1.00  msec 

The  time  response  is  formed  using  only  the  lowest  3 natural  frequencies 
and  corresponding  normal  jnodes. 

TABLE  2:  Numerical  model  simulation  of  a UHF  Polemast  Antenna.  All  symbols 
are  defined  in  the  text.  The  structural  data  were  obtained  from 
drawings  supplied  by  DMCS-6  (modifications  as  noted  in  [13]).  The 
air  blast  predictions  were  provided  by  the  United  States  Defense 
Nuclear  Agency.. 
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Mode 

Whip  Frequencies  (cps) 

Polemast  Frequencies  (cps) 

Theoretical 

Experimental 

Theoretical 

Experimental 

1 

1.47 

1.27 

4.62 

4.00 

2 

4.09 

4.20 

25.5 

24. l’ 

3 

9.55 

9.50 

72.4 

— 

This  value  represents  an  average  of  indistinct  frequencies  which 
appear  in  a band  over  the  range  19.7  to  32.1  cps. 


TABLE  3:  A comparison  of  theoretical  (numerical  simulation)  and 

experimental  (Twang  Test)  natural  frequencies  for  a 35  ft 
fibreglass  Whip  Antenna  and  a 23  ft  UHF  Polemast  Antenna. 
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FIGURE  1: 


Photograph  of  the  35  ft  fibreglass  Whip  Antenna  installed  at 
the  nominal  7.0  psi  location  iTi  Event  Dice  Throw. 
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FIGURE  3:  Comparison  of  the  Whip  Antenna  bending  strain  predictions 
(dashed  lines)  against  the  measured  strains  (solid  lines). 

The  five  gauge  pairs  are  respectively  located  at  distances 
3.5,  10. S,  17.0,  18.4  and  24.0  ft  from  the  base  of  the  antenna. 
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FIGURE  4:  Comparison  of  the  Polemast  Antenna  bending  strain  predictions 
(dashed  lines)  against  the  measured  strains  (solid  lines). 

The  five  gauge  pairs  are  respectively  located  at  distances 
0.?5,  1.5,  3.Z5,  4.5  and  12.0  ft  from  the  base  of  the  antenna. 
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NC]MINAL  VS-  EXPERIMENTAL  OVERPRESSURE 

FiajRE  5:  Comparison  of  the  Friedlander  wave  which  corresponds  to  the  pre-trial  DNA 
overpressure  predictions  against  the  Friecllander  wave  which  corresponds 
to  average  experimental  overpressure  measurements. 
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APPENDIX  A 

COMPUTER  PROGRAM  FOR  THE  NUMERICAL  SIMULATION  MODEL 

Based  on  the  numerical  procedures  described  in  Sections  2,  3 and  4, 
a computer  program  war.  prepared  to  determine  the  natural  frequencies,  normal 
modes  and  air  blast  transient  response  of  a tapered  cantilever  beam.  A 
listing  of  the  computer  program  is  presented  in  this  appendix,  together  with 
a description  of  the  initialization  procedures. 

(a)  Free  Vibration 

The  computer  program  for  free  vibration  consists  of  a main  program 
and  two  subroutines,  referred  to  as  "ESOR"  and  "ANORM".  The  major  portion  of 
the  computation  for  the  natural  frequencies  and  nonnal  modes  is  carried  out 
in  subroutine  ESOR.  On  completion  of  the  calculation,  the  program  displays 
the  relevant  frequencies  and  modes,  and  stores  the  information  in  a disk  file 
for  future  use  in  the  forced  vibration  calculations.  The  card  input  data  for 
the  free  vibration  calculation  is  outlined  below,  and  a listing  of  the 
computer  program  is  provided  in  Figure  A1 . 
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Program 

Symbol 

Text 

Symbol 

Units 

Field 

Type 

Description 

CARD  1;  Iteration  control  parameters 

N 

N 

— 

15 

number  of  elements  in  the  beam. 

M 

15 

number  of  frequencies/modes  to 
be  calculated,  ascending  order 
(3  is  satisfactory). 

FSTOP 

F10.5 

SOR  and  Rayleigh  quotient  con- 
vergence criteria  (.001  is 
satisfactory. 

NSTOP 

- 

— 

15 

iteration  loop  termination 
parameter  (150  is  satisfactory). 

IPRNT 

- 

— 

15 

printout  option  (1=full, 

0=partial ) . 

ALEN 

L 

-- 

F10.5 

length  of  the  beam. 

NN 

15 

number  of  data  cards  (see  Card 

2)  used  to  describe  the  mass 
distribution  along  the  antenna. 

EMOD 

E 

psi 

FI  0.5 

Young's  modulus. 

RHO 

P 

slugs/in^ 

F10.5 

density. 

ANU 

- 

-- 

FID.  5 

son  acceleration  factor  (1.0  is 
sati sfactory) . 

IBC 

I 

1 

15 

boundary  condition  code: 
0=clamped-free,  1 or  larger  = 
pin-pin-free  with  pin  at 

CARD  2;  Description  of  the  mass  distribution  along  the  beam;  NN  cards  of 
this  type. 


Wl.XX 

X 

ft 

Flo. 2 

position  from  base  of  the  beam. 

mSM 

in 

FI  0.2 

outside  diameter. 

m 

ID 

in 

F10.2 

inside  diameter. 
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FREE  VIBRATION  (CONT’D) 
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29. 


(b)  Fot'ced  Vibration 

The  computer  program  for  forced  vibration  consists  of  a main 
program  and  three  subroutines,  referred  to  as  "DRAGC"*  "YDISP"  and  "VISCO". 
The  major  portion  of  the  computation  for  the  forced  vibration  response  is 
carried  out  in  the  main  program.  Subroutines  DRAGC,  YDISP  and  VISCO 
respectively  deal  with  computation  of  the  air  blast  loading  on  circular 
cylinders,  display  of  the  stresses  which  correspond  to  a specified  beam 
displacement  profile,  and  calculation  of  the  kinematic  viscosity  for  air. 
Input  data  for  this  program  consists  of  the  natural  frequencies  and  normal 
modos  computed  previously  (read  from  a disk  file)  and  data  cards  identifying 
the  air  blast  loading  and  assorted  parameters  required  in  the  time  response 
calculation.  The  card  input  data  for  the  forced  vibration  calculation  is 
outlined  below,  and  a listing  of  the  computer  program  is  provided  in  Figure 
A2. 


Program 

Symbol 

Text 

Symbol 

Units 

Field 

Type 

Description 

Card  1 ; 

Data  describing  the 

normal  mode  plots. 

INPLT 

- 

-- 

15 

plotting  option  (0=5uppress  the  plot). 

ICODE 

- 

15 

printout  option  in  the  plot  subroutine. 

XLEN 

- 

in 

F10.3 

length  of  the  abscissa. 

YUEN 

- 

in 

Flo. 3 

length  of  the  ordinate. 

YDELL 

~ 

— 

Flo. 3 

scale  of  the  ordinate. 

UNCLASSIFIED 


UNCLASSIFIED 


30. 


Text 

Symbol 

Units 

Field 

Type 

Description 

Card  2:  Data  describing  the  air  blast  loading  on  the  antenna.  | 

PO 

Po 

psi 

FlO.l 

peak  overpressure. 

TPLUS 

^d 

sec 

FlO.l 

positive  phase  duration. 

PA 

psi 

FlO.l 

atmospheric  pressure. 

TA 

Ta 

op 

FlO.l 

atmospheric  temperature. 

D 

in 

FlO.l 

representative  antenna  diameter  (for 
display  of  drag  pressures  only). 

AK 

K 

FlO.l 

Fried iander  decay  constant. 

TIME 

- 

sec 

FlO.l 

time  step  (for  display  of  drag  pressures 
only;  .010  is  satisfactory). 

LOOPS 

L 

no 

limit  on  the  number  of  display  loops 
(display  purposes  only;  999  is 
satisfactory) . 

Card  3:  Number  of  cards  needed  to  specify  the  projected  (notmal  to  blast 

direction) 

cross-sectional  area  profile  of  the  antenna.  1 

NN 

15 

if  NN-1 j tnen  tne  OD  used  in  the  Free 
Vibration  problem  also  specifies  the 
projected  antenna  area  profile;  if 

NN>1 , then  there  are  NN  cards  to  follow 
to  specify  this  profile  (see  Card  4). 

Card  4:  Description  of  the  projected  (normal  to  blast  direction)  cross-sectionaili 

-.-.iCaa.ECQfi' 



W1  .XX 

X 

ft 

Flo. 2 

position  from  base  of  the  beam. 

W3,DIA 

OD 

in 

Flo. 2 

outside  diameter. 

1 

|Card  5:  Control  par< 

imeters  foi 

" the  tran 

>ient  response  calculation. 

MODES 

no 

number  of  normal  modes  used  in  construc- 
ting the  time  response  (ascending  order 
in  frequency). 

DELT 

At 

sec 

no.  5 

time  step  used  in  the  numerical 
integration. 

TLIM 

- 

sec 

no 

maximum  time  to  which  the  time  response 
is  carried. 

MDISP 

F1C.5 

the  number  of  time  steps  between  beam 
response  displays. 

UNCLASSIFIED 
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